Integer quantum Hall transition in the presence of a long-range-correlated quenched 

disorder 
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We theoretically study the effect of long-ranged inhomogeneities on the critical properties of the 
integer quantum Hall transition. For this purpose we employ the real-space renormalization-group 
(RG) approach to the network model of the transition. We start by testing the accuracy of the RG 
approach in the absence of inhomogeneities, and infer the correlation length exponent v = 2.39 from 
a broad conductance distribution. We then incorporate macroscopic inhomogeneities into the RG 
procedure. Inhomogeneities are modeled by a smooth random potential with a correlator which falls 
off with distance as a power law, r~" . Similar to the classical percolation, we observe an enhancement 
of u with decreasing a. Although the attainable system sizes are large, they do not allow one to 
unambiguously identify a cusp in the v{a) dependence at = Ijv, as might be expected from 
the extended Harris criterion. We argue that the fundamental obstacle for the numerical detection 
of a cusp in the quantum percolation is the implicit randomness in the Aharonov-Bohm phases of 
the wave functions. This randomness emulates the presence of a short-range disorder alongside the 
smooth potential. 



FAGS numbers; 73.40.Hm, 61.43.-j, 64.60.Ak 
I. INTRODUCTION 

The critical behavior of electron wave functions in the 
vicinity of the integer quantum Hall (QH) transition is 
now well understood^ That is, the locaUzation length di- 
verges as where e is the deviation from the critical 
energy. The most accurate value of the exponent v exj-i 
tracted from numerical simulations is = 
On the experimental side, the study of the critical be- 
havior of the resistance in the transition region at strong 
magnetic field B has a long history which can be conven- 
tionally divided into three periods. _ _ 

(a) The first experimental worksOij reported a nar- 
rowing of the transition peak with temperature T, as 
with K ~ 0.4. The spread in the actual value of k oc 1/;/ 
measured in different experiments was attributed to the 
difference in the type of disorder in the samples of Refs. 

and ||. Another experimental method to explore the 
critical behavior was employed in Refs. |^ and |^ where 
K was deduced from the sample size dependence of the 
width, AB, of the transition region. The value of n ob- 
tained by this technique appeared to be consistent with 
temperature measurements of Ref. ^, in the sense that n 
was found to be sample dependent. On the other hand, it 
was argued in Ref. ||that the lack of universality in Refs. 

has its origin in the long-ranged character of the dis- 
order in GaAs-based heterostructures studied in these 
works. This is because for a smooth disorder the energy 
interval within which the electron transport is dominated 
by localization effects is relatively narrow.B The measure- 
ments in Refs. ^ and ^ suggesting the universality of k 
were carried out on In2:Gai_2;As/InP heterostmptures in 
which disorder is believed to be short-ranged .E2l Despite 



the disagreement about universality of the exponent n, 
the fact that the narrowing of the plateau transition oc- 
curs as T'^ was not questioned in Refs. 

(b) The absence of scaling was reported first for 
the QH-insulator tpansitiontll and then for the plateau- 
plateau transition.E3 In the latter paper the conclu- 
sion about the absence of scaling was drawn from 
the analysis of the frequency dependence of AB in 
GaAs/Alj^Gai_j,As heterostructures (in contrast to the 
similar analysis in Ref. That is, the authors of Refs. 
pd| and |l^ concluded, that the width of the transition 
region saturates as T — > 0. A possible explanation of 
this behaviortSEj is based on the scenario of tunneling 
between electron puddles with a size larger than the de- 
phasing length. The microscopic origin ofpt^psp puddles 
was attributed to sample inhomogeneitie 

(c) Very recent experimental resultal; 
plateau-insulator as well as plateau-plateau QH transi- 
tions carried out on the same Inj^Gai-^^As/InP sample as 
in Ref. || suggested that the narrowing of the transition 
width with temperature follows a power-law dependence 
AB (X T'^ with K « 0.4. Even when the authors of Ref. 
|T^ analyzed their data using the procedure of Ref. |l^, 
i.e., by plotting the logarithm of the longitudinal resis- 
tance versus AB, they obtained straight lines with slopes 
proportional to with k' w 0.55. They attributed the 
difference between k and k' to the marginal dependence 
of the critical resistance on T. It was also speculated 
in Ref. |l^ that this dependence most likely results from 
macroscopir inhomogeneities in the sample. In the latest 
papersE2rE3 the frequency dependence of the QH transi- 
tion width was studied. Ihe-iresults did not support the 
saturation of the width,MO but rather confirmed the 
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scaling hypothesis. 

Summarizing, it is now conclusively established that 
insulator-plateau and plateau-plateau transitions exhibit 
the same critical behavior. It is also recognized that 
macroscopic inhomogepeities can either mask the scaling 
or alter the value of kO 

On the theoretical side, in all previous considera- 
tions inhomogeneities were incorporated into thft-Uieops 
through a spatial variation of the local resistivityliiliiSa 
In other words, all existing theories are either "homo- 
geneous quantum coherent" or " inhomogeneous incoher- 
ent . Meanwhile, there is another scenario which has 
never been explored. Close to the transition the quantum 
localization length ^ becomes sufficiently large. Then the 
long-ranged disorder can affect the character of the diwjs-. 
gence of ^. At this point we recall the classical limit ,cJ'E3 
in which the long-ranged disorder does affect the value 
of the critical exponent in the percolation problem. Ob- 
viously, when the disorder is long-ranged but has a finite 
correlation radius, one should not expect any changes in 
the critical behavior. The principle finding of Refs. |2j 
and ^ is that the critical exponent can change when 
the correlator of the disorder {V{r)V{r')) falls off with 
distance as a power law, i.e., oc |r — r'|~", (quenched 
disorder). According to Refs. |2j and ^ the critical ex- 
ponent of the classical percolation = 4/3 crosses over 
to = 2/a for a < 3/2, i.e., when the decay of the corre- 
lator is slow enough. In the present paper we study the 
effects of quenched disorder on quantum percolation. The 
latter is known to describe the localization-delocalization 
transition for a two-dimensional (2D) electron in a strong 
magnetic field. As a model of quantum percolation we 
employ the Chalker-Coddington (CC) modelcj which is 
one of the mainJ' tools" for the quantitative study of the 
QH transition.EZrH The CC model is a strong magnetic 
field (chiral) limit of.^ general network model, first in- 
troduced by ShapircEj and later utilized for the study 
of localization-delocalization transitions within different 
universality classes .e3^c£| In addition to describing the 
QH transition, the CC model applies to a much broader 
class of critical phenomena since the correspondence be- 
tween the CC model and_themiodynamic, field-theory 
and Dirac-fermions modelsOLj was demonstrated. 

In order to study the role of the quenched disorder on 
the localization-delocalization transition we treat the CC 
model within the real-space renormalization group (RG) 
approach. First, in Sec. || we check the accuracy of the 
RG approach, and show that it provides a remarkably 
accurate description of the QH transition. In Sec. 
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we extend the RG approach to incorporate the quenched 



disorder. Concluding remarks are presented in Sec. IV 



II. TEST OF THE RG APPROACH TO THE CC 
MODEL 



A. Description of the RG approach 

As mentioned in Sec. |, thajCC model is a chiral limit of 
the general network model.Lj However, it was originally 
derived from a microscopic picture of electron motion in 
a strong magnetic field and a smooth potential.E3 Within 
this picture, the links can be identified with semiclassical 
trajectories of the guiding centers of the cyclotron orbit, 
while the nodes correspond to the saddle points (SP's) at 
which different trajectories come closer than the Larmour 
radius. For simplicity, the nodes were placed on a square 
lattice. We will use an equivalent, but slightly different 
graphical representation of the CC network shown in Fig. 
|l. In this representation the centers of the trajectories of 
the guiding centers are placed on a square lattice and play 
the role of nodes, whereas the SP's should be identified 
with links. r-ir-i 

We now apgly; a real-space RG approacfSEj to the 
CC network.E3E3 The RG approach is based on the as- 
sumption that a certain part of the network containing 
several SP's, the RG unit, represents the entire network. 
This unit is then replaced through the RG transforma- 
tion by a single super-SP with an S matrix determined by 
the S matrices of the constituting SP's. The network of 
super-SP 's is then treated in the same way as the original 
network. Successive repetition of the RG transformation 
yields the information about the S matrix of very large 
samples, since, after each RG step, the effective sample 
size grows by a certain factor determined by the geome- 
try of the original RG unit. Obviously, a single RG unit 
is a rather crude approximation of the network. There- 
fore, prior to applying the RG approach to the study of 
the quantum Hall transition in the presence of quenched 
disorder, we first check the accuracy of this approach for 
the conventional, uncorrelated case, where the compar- 
ison with the results of direct numerical simulations is 
possible. 

The RG unit we use is extracted from a CC network on 
a regular 2D square lattice as shown in Fig. |l|. A super- 
SP consists of five original SP's by analogy to the[R|(i-i|init 
employed for the 2D bond percolation problem.Ea^Ell As 
in any RG scheme, the unit shown in Fig. |] leaves out a 
number of bonds of the original lattice. Nevertheless, it 
is well known that the application of this-scheme to the 
classical case yields very accurate resultsO 

Between the SP's an electron travels along equipo- 
tential lines, and accumulates a certain Aharonov-Bohm 
phase. Different phases are uncorrelated, which reflects 
the randomness of the original potential landscape. Each 
SP can be described by two equations relating the wave- 
function amplitudes in incoming and outgoing channels. 
This results in a system of ten linear equations, the so- 
lution of which yields the following expression for the 
transmission coefficient of the super-SP (Ref. pq): 
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tMr2r3rie'^' ~ 1) + ^2^4e^^^^+'*'"H^i^3r5e-"*'^ - 1) + ^3(^2^56^^^ + hUe"^^) 



Here ti and = (1— )^/^ are, respectively, the transmis- 
sion and reflection coefficients of the constituting SP's; 
$j are the phases accumulated along the closed loops (see 
Fig. I). Equation (0) is the RG transformation, which 
allows one to generate (after averaging over <i>j ) the dis- 
tribution P{t') of the transmission coefficients of super- 
SP's using the distribution P{t) of the transmission coef- 
ficients of the original SP's. Since the transmission coef- 
ficients of the original SP's depend on the electron energy 
£, the fact that delocalization occurs at e = implies that 
a certain distribution, Pc{t) — with Pc{t^) being symmet- 
ric with respect to = — is the fixed point (FP) of 
the RG transformation Eq. ([!]). The distribution Pc{G) 
of the dimensionless conductance G can be obtained from 
the relation G = i^, so that Pc{G) = Pc{t)/2t. 

B. Critical exponent within the RG approach 

Since the dimensionless SP height Zi and the trans- 
mission coefficient at e = are related as ti = 
(e^' + transformation (|^) determines the height 

of a super-SP through the heights of the five constitut- 
ing SP's. Correspondingly, the distribution P{G) de- 
termines the distribution Q{z) of the SP heights via 
Q{z) = P{G){dG/dz) = icosh-2(z/2)P[(e" + 1)-^]. 
In fact, Q{z) is not a characteristic of the actual SP's, 
but rather, as we will see below, represents a convenient 
parametrization of the conductance distribution. 

The language of the SP heights provides a natural 
way to extract the critical exponent v. Suppose that 
the RG procedure starts with an initial distribution, 
Qq{z) — Qc{z — zq), that is shifted from the critical dis- 
tribution, Qdz), by a small zq oc e. Since zq <C 1, the 
first RG step would yield Qc{z — tzq) with some num- 
ber T independent of zq. After n steps the center of the 
distribution will be shifted by Zmax.n = t"'zq, while the 
sample size will be magnified by 2". At the nth step 
corresponding to Zmax,n ~ 1 a typical SP is no longer 
transmittable. Then the localization length ^ should be 
identified with 2" oc Zq'^ oc e^'^, with v = In 2/ In t. 
When the RG procedure is carried out numerically, one 
should check that Zq is small enough so that Zmax,n oc zq 
for large enough n. Consequently, the working formula 
for the critical exponent can be presented as 



which should be independent of n for large n. 



C. Numerical results 

In order to find the FP conductance distribution 
Pc{G), we start from a certain initial distribution of 
transmission coefficients, Po{t) (see below). The distri- 
bution is discretized in at least 1000 bins, such that the 
bin width is typically 0.001 for the interval t € [0,1]. 
From ^0(^)7 we obtain ti, i — 1,...,5, and substitute 
them into the RG transformation [Eq. (PJ. The phases 
<I>j, j ~ 1, ... ,4 are chosen randomly from the interval 
e [0, 27r]. In this way we calculate at least 10^ super- 
transmission coefficients t'. The obtained histogram 
Pi(t') is then smoothed using a Savitzky-Golay filterH in 
order to decrease statistical fiuctuations. At the next step 
we repeat the procedure using Pi as an initial distribu- 
tion. Wc assume that the iteration process has converged 
when the mean-square deviation J dt[Pn{t) — Pn-i{t)]'^ 
of the distribution P„ and its predecessor P„_i deviate 
by less than 10"''. 

We are now able to study samples with short-ranged 
disorder. The actual initial distributions, Po(i), were 
chosen in such a way that corresponding conductance 
distributions, Po(G), were either uniform or parabolic, 
or identical to the FP distribution found semianalyti- 
cally in Ref. All these distributions are symmetric 
with respect to G = 0.5. We find that, regardless of the 
choice of the initial distribution, after 5-10 steps the RG 
procedure converges to the same FP distribution which 
remains unchanged for another 4-6 RG steps. Small de- 
viations from symmetry about G = 0.5 finally accumu- 
late due to numerical instabilities in the RG procedure, 
so that typically after 15-20 iterations the distribution 
becomes unstable and flows towards one of the classical 
FP's P(G) = (5(G) or P(G) = S{G - 1). We note that 
the FP distribution can be stabilized by forcing Pn{G) 
to be symmetric with respect to G = 0.5 in the course of 
the RG procedure. 

Figure |2| illustrates the RG evolution of P(G) and 
Q{z). In order to reduce statistical fluctuations we aver- 
age the FP distributions obtained from different Po(G)'s. 
The FP distribution Pc(G) exhibits a flat minimum 
around G = 0.5, and sharp peaks close to G = and 
G = 1. It is symmetric with respect to G « 0.5 with 
(G) = 0.498 ± 0.004, where the error is the standard er- 
ror of the mean of the obtained FP distribution. The FP 
distribution Qc{z) is close to Gaussian. 

We now turn to the critical exponent u. As a result of 
the general instability of the FP distribution, an initial 
shift of Qc{z) by a value zq results in the further drift 
of the maximum position, 0max,n, away from z = af- 
ter each RG step. As expected, Zmax,™ depends linearly 
on Zq. This dependence is shown in Fig. ^ (inset) for 
different n from 1 to 8. The critical exponent is then 
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calculated from the slope according to Eq. (||). Figure ^ 
illustrates how the critical exponent converges with n to 
the value 2.39 ± 0.01. The error corresponds to a confi- 
dence interval of 95% as obtained from the fit to a linear 
behavior. 

Due to the high accuracy of our calculation of Pc{G), 
we were able to reliably determine many central moments 
((G - (G))") of the FP distribution Pc(G). These mo- 
ments are plotted in Fig. 0. 



D. Comparison with previous simulations 

By dividing the CC network into units, the RG ap- 
proach completely disregards the interference of the 
wave-function amplitudes between different units at each 
RG step. For this reason it is not clear to what extent this 
approach captures the main features and reproduces the 
quantitative predictions at the QH transition. Therefore, 
a comparison of the RG results with the results of direct 
simulations of the CC model is crucial. These direct sim- 
ulations are usualljy carried out by enApJoying either the 
quasi- ID versioiB or the 2D versionEJ of the transfer- 
matrix method. The results of application of the version 
of Ref. |6^ to the CC model are reported in Refs. ^ and 
^t\ . In Refs. |65|-|67|the version of Ref. |64|was utilized. For 
the critical exponent the values v — 2.5 ± 0.5 (Ref. |2^ ) 
and later v = 2.4±0.2 (Ref. 27) were obtained. Note that 
our result is in excellent agreement with these values, and 
is also consistent with the most precise = 2.35 ± 0.03.0 
This already indicates a remarkable accuracy of the RG 
approach. In Refs. |65|"|67| the critical distribution Pc{G) 
of the conductance was studied. Pc{G) was found to 
be broad, which is in accordance with Fig. ||. However, 
a more detailed compafdsmi is impossible, since the re- 
sults of the simulationsE3~E3 do not obey the electron-hole 
symmetry condition Pc{G) — Pc{l — G). On the other 
hand, within the RG approach, the latter condition is sat- 
isfied automatically. Nevertheless, we can compare the 
moments of Pc{G) to those calculated in Refs. ^ and |6^. 
They are presented in Fig. |[ In Ref. ^ only the stan- 
dard deviation ((G^) — (G)^) "^^^ « 0.3 was computed. 
Our result is 0.316. In Ref. p5^the moments were fitted 
by two analytical functions, which are also shown in Fig. 
^. They agree with our calculations up to the sixth mo- 
ment. Here we point out that the moments obtained in 
Ref. ^ can hardly be distinguished from the moments of 
a uniform distribution. This reflects the fact that Pc{G) 
is practically flat except for the peaks close to G = and 
G = 1. 

In Refs. ^ and|69|Pc(G) was studied by methods which 
are not based on the CC model. Both works reported a 
broad distribution Pc(G). In Ref. ^Pc(G) was found to 
be almost flat. The major difference between Ref. |6^ and 
Fig. ||is the behavior of Pc{G) near the points G = and 
1. That is, P(G) drops in Ref. |6| to zero at the ends, 
while Fig. exhibits maxima. In Ref. 69 the behavior 



of Pc (G) is qualitatively similar to Fig. g, with maxima 
at G = and 1. However, the statistics in Ref. ^ are 
rather poor, which again rules out the possibility of a 
more detailed comparison with our results. 

Finally, we p oint out that our results agree completely 
with Refs. [70|-|72| where a similar RG treatment of the 
CC model was carried out. Our numerical data have a 
higher resolution, and show significantly less statistical 
noise. This is because we took advantage of faster com- 
putation by using the analytical solution of the RG [Eq. 
(|l|)].Ea Also, note that, in Refs. ^ and |7| the critical ex- 
ponent v was calculated using a procedure different from 
that described in Sec. II B, Nevertheless, the values of 



ly determined by both methods are close. We emphasize 
that a systematic improvement of the RG procedure, i.e., 
by inclusion of more than five SP's into the basic RG unit 
as reported in Refs. [70[-[72|, leads to similar results. The 
critical distributiop_p£ the conductance was also stud- 
ied experimentallylljo in mesoscopic QH samples. Al- 
though an almost uniform conductance distribution con- 
sistent with theoretical predictions was found in Rje£ |7^, 
further detailed analysis of the mesoscopic patternO has 
revealed the crucial role of the charging effects, which 
were neglected in all theoretical studies. 

We conclude that the test of the RG approach against 
other simulations proves that this approach provides a 
very accurate quantitative description of the QH transi- 
tion. It can now be used to study the influence of long- 
ranged correlations. 



III. MACROSCOPIC INHOMOGENEITIES 

A. General considerations 

A natural way to incorporate a quenched disorder into 
the CC model is to ascribe a certain random shift, zq, to 
each SP height, and to assume that the shifts at different 
SP positions, r and r', are correlated as 



{zQir)zQ{r')) oc 



(3) 



with a > 0. After this, the conventional transfer-matrix 
methods of Refs. ^ and ^ could be employed for nu- 
merically precise determination of (G), the distribution 
Pc(G), its moments, and most importantly, the critical 
exponent, ly. However, the transfer-matrix approach for 
a 2D sample is usually limited to fairly small sizes (e.g., 
up to 128 in Ref. ^) due to the numerical complexity 
of the calculation. Therefore, the spatial decay of the 
power-law correlation by, say, more than an order in mag- 
nitude is hard to investigate for small .q;. In principle the 
quasi- ID transfer-matrix methocOO can easily handle 
such decays at least in the longitudinal direction, where 
typically a few million lattice sites are considered iter- 
atively. A major drawback, however, is the numerical 
generation of power-law correlated randomness since no 
iterative algorithm is knownJl3ll3 This necessitates the 
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complete storage p£ different samples of correlated SP 
height landscapesJIZi and the advantage of the iterative 
transfer-matrix approach is lost|— .Furthermore, in order 
to deduce the critical exponent ,EZl one needs to perform 
finite-size scaling£3 with transverse sizes that should also 
be large enough to capture the main effect of the power- 
law disorder in the transverse direction. Consequently, 
even for a single transfer-matrix sample, the memory re- 
quirements add up to gigabytes. 

On the other hand, the RG approach is perfectly suited 
to study the role of the quenched disorder. First, after 
each step of the RG procedure, the effective system size 
doubles. At the same time, the magnitude of the smooth 
potential, corresponding to the spatial scale r, falls off 
with r as r""/^. As a result, the modification of the 
RG procedure due to the presence of the quenched dis- 
order reduces to a rescaling of the disorder magnitude 
by a constant factor 2~"/^ at each RG step. Second, 
the RG approach operates with the conductance distri- 
bution, Pn{G), which carries information about all the 
realizations of the quenched disorder within a sample 
of a size 2^|— .This is in contrast to the transfer-matrix 
approach, E^lEil within which a small increase of the sys- 
tem size requires the averaging over the quenched disor- 
der realizations to be conducted again. 

In order to find out whether or not the critical behavior 
is affected by the quenched disorder, the following argu- 
ment was put forward in Ref. In the absence of the 
quenched disorder, the correlation length, for a- given 
energy zq in the vicinity of the transition is proportional 



to z, 



Q 



Now consider a sample with an area A = ^q. 
The variance of the quenched disorder within the sample 
is given by 



^0 



= i^{lAd'rzQ{r)J^d\'zQ{r')) 
= wlAd'rhd^r'{zQir)zQ{r')} 



(4) 



where the last relation follows from Eq. (||). 

The critical behavior remains unaffected by the 
quenched disorder if the condition Aq/zq — > as zq — > 
is met. Using Eq. the ratio Aq/zq can be presented 




ln(zQ'^), a — 2 
a<2 



(5) 



We thus conclude that quenched disorder is irrelevant 
when 



v > I, for a > 2, 
ai' > 2, for a < 2 



(6) 



The first-condition corresponds to the original Harris 
criterion^ for uncorrelated disorder, while the second 
condition is the extended Harris criterion.cil It yields the 
critical value of the exponent a, i.e., ac = 2/ v. 



The above consideration suggests the following algo- 
rithm. For the homogeneous case all SP's constituting 
the new super-SP are assumed to be identical, which 
means that the distribution, of heights, Qniz), for all of 
them is the same. For the correlated case these SP's are 
no longer identical, but rather their heights are randomly 
shifted by the long-ranged potential. In order to incorpo- 
rate this potential into the RG scheme, Qn{zi) should be 

replaced by Qn{zi — A|"^) for each SP, i, where A^"-* is 
the random shift. Then the power-law correlation of the 
quenched disorder enters into the RG procedure through 
the distribution of A^p\ That is, for each n the distri- 
bution is Gaussian with the width /3(2")~"/^. For large 
enough n the critical behavior should not depend on the 
magnitude /3, but on the power, a, only. 



B. Numerical results 

Here we report the results of the application of the al- 
gorithm outlined in the previous section. First, we find 
that for all values of a > in correlator the FP dis- 
tribution is identical to the uncorrelated case within the 
accuracy of our calculation. In particular, (G) = 0.498 is 
unchanged. However, the convergence to the FP is nu- 
merically less stable than for uncorrelated disorder due to 
the correlation-induced broadening of Qniz) during each 
iteration step. In order to compute the critical exponent 
i/(a) we start the RG procedure from Qo{z — zo), as in the 
uncorrelated case, but, in addition, we incorporate the 
random shifts caused by the quenched disorder in gener- 
ating the distribution of z at each RG step. The results 
shown in Fig. [s] illustrate that for increasing long-ranged 
character of the correlation (decreasing a) the conver- 
gence to a limiting value of v slows down drastically. 
Even after eight RG steps (i.e., a magnification of the 
system size by a factor of 256), the value of v with long- 
ranged correlations still differs appreciably from ly — 2.39 
obtained for the uncorrelated case. The RG procedure 
becomes unstable after eight iterations, i.e., Zmax.g can no 
longer be obtained reliably from Qq{z). Unfortunately, 
for small a the convergence is too slow to yield the limit- 
ing value of v after eight steps only. For this reason, we 
are, strictly speaking, unable to unambiguously answer 
the question whether sufficiently long-ranged correlations 
result in an a-dependent critical exponent J^(a), or the 
value of v eventually approaches the uncorrelated value 
of 2.39. Nevertheless, the results in Fig. |^ indicate that 
the effective critical exponent exhibits a sensitivity to the 
long-ranged correlations even after a large magnification 
by 256x256. Therefore, in realistic samples of finite sizes, 
macroscopic inhomogeneities are able to affect the results 
of scaling studies. Note further that as shown in the in- 
set of Fig. ^ there is no simple scaling of values when 
plotted as function of an rcnormalized system size 2""/^. 
We emphasize that i/(q;) obtained after eight RG steps 
always exceeds the uncorrelated value. Thus, our results 
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indicate that macroscopic inhomogeneities must lead to 
smaller values oi k (xl/v. Experimentally, the value of k 
smaller than 0.42 was reported in a number of early (see, 
e.g., Ref. and references therein) as well as recentL3 
works. This fact was accounted for by different reasons 
(such as temperature dependence of the phase breaking 
time, incomplete spin resolution, valley degeneracy in Si- 
based metal-oxide-semiconductor field-effect transistors, 
and inhomogeneity of the carrier concentration in GaAs- 
based structures with a wide spacer). Briefly, the spread 
of the K, values was attributed to the fact that the tem- 
peratures were not low enough to assess the truly critical 
regime. The possibility of having k < 0.4 due to the 
correlation-induced dependence of the effective v on the 
phase-breaking length or, ultimately, on the sample size, 
as in Fig.Ja, was never considered. 

Figure |[shows the values of v obtained after the eighth 
RG step as a function of the correlation exponent a for 
different dimcnsionless strengths /3 of the quenched dis- 
order. It is seen that in the domain of a, where the values 
of v differ noticeably from v = 2.39, their dependence on 
(3 is strong. According to the extended Harris criterion, 
v{q) is expected to exhibit a cusp at the /3-independent 
value a = ac = 2/2.39 w 0.84. From our results in 
Fig. H, two basic observations can be made. For a small 
enough magnitude of the long-range disorder, we see a 
smooth enhancement of v{a) with decreasing a without 
a cusp. Although such a behavior is due to the relatively 
small number of RG steps, the data might be relevant 
for realistic samples which have a finite size and a finite 
phase-breaking length governed by temperature. At the 
largest /3 = 4, the cusp eventually shows up but the nu- 
merics becomes progressively ambiguous, forbidding us 
from going to even larger /3. The origin for this strong 
(3 dependence of our results is a profound difference be- 
tween the classical and quantum percolation problems. 
This difference is discussed in the Sec. 



IV. DISCUSSION 

A. Classical case 

In the classical limit, the motion of an a electron in 
a strong magnetic field and a smooth potential reduces 
to the drift of the Larmour circle along the equipotential 
lines. Correspondingly, the description of the delocaliza- 
tion transition reduces to the classical percolation prob- 
lem. As mentioned above, for classical percolation the 
quenched disorder is expected to cause a crossover in the 
exponent Vc, describing the size of a critical equipotential 
from Vc = 4/3 to Vc = 2/a for a < 3/2. This predictionEa 
was made on the basis-pf Eq. (^). It was later tested by 
numerical simulationsL3 which utilized the Fourier filter- 
ing method to generate a long-range-correlated random 
potential. The exponent vdo/) waS|-S|tudied using the 
same classical real-space RG methodEJ that we utilized 



above. The values of Vc inferred for a > 1 were consis- 
tently lowet-^han 2/a. For example, Vc — ^ was found 
for a — 0.4,c3 whereas Vc — 3.4 ± 0.3 was observed. 

The classical version of the delocalization transition is 
instructive, since it allows one to trace how the critical 
equipotentials grow in size upon approaching the percola- 
tion threshold, and how the quenched disorder affects this 
growth. Roughly speaking, in the absence of long-ranged 
correlations, the growth of the equipotential size is due 
to the attachment of smaller equipotentials to the criti- 
cal ones. As a result, the shape of a critical equipotential 
is dendritelike. As the threshold is approached, differ- 
ent critical equipotentials become connected through the 
narrow "arms" of the dendrite. Long-ranged correlations 
change this scenario drastically. As could be expected 
intuitively, and as follows from the simulations criti- 
cal equipotentials become more compact due to correl^ 
tions. In fact, for a < 0.25, the "arms" play no rolc,LJ 
i.e., the morphology of a critical equipotential becomes 
identical to its "backbone". As a result, the formation of 
the infinite equipotential at the threshold occurs through 
a sequence of "broad" merges of compact critical equipo- 
tentials. The correlation-induced enhancement of Vc in- 
dicates that due to these merges the size of the critical 
equipotential in the close vicinity of the threshold grows 
faster than in the uncorrelated case. Since our simula- 
tions also demonstrate the enhancement of the critical 
exponent due to the correlations, the main result of the 
present paper can be formulated as follows: quenched 
disorder affects classical and quantum percolation in a 
similar fashion. 



B. Quantum case 

Here we note that there is a crucial distinction between 
the classical case and the quantum regime of the electron 
motion considered in the present paper. Indeed, within 
the classical picture, correlated disorder implies that the 
motion of the guiding centers of the Larmour orbits in 
two neighboring regions is completely identical. In our 
study, we have incorporated the correlation of heights 
of the saddle points into the RG scheme. At the same 
time we have assumed that the Aharonov-Bohm phases 
acquired by an electron upon traversing the neighbor- 
ing loops are completely uncorrelated. This assumption 
implies that, in addition to the long-ranged potential, 
a certain short-ranged disorder causing a spread in the 
perimeters of neighboring loops of the order of the mag- 
netic length is present in the sample. The consequence of 
this short-range disorder is the sensitivity of our results 
to the value of /3 which parametrizes the magnitude of 
the correlated potential. The presence of this short-range 
disorder affecting exclusively the Aharonov-Bohm phases 
significantly complicates the observation of a cusp in the 
j^(a) dependence at a « 0.84, as might be expected from 
the extended Harris criterion. 
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Let us elaborate on these complications. A general 
form of the correlator for long-range disorder is 



{V{r)V{r')) = V^F 



|r-r'| 
I 



(7) 



where I is the microscopic length; -F(O) — 1, F[x) oc x^" 
for a; S> 1. Now suppose that the correlator contains 
an additional short-range term W^G(|r — r'|/^), with 
G{0) = 1 and G{x) falling off much more rapidly than 
F{x) for a; » 1. The extended Harris criterion implies 
that this term will not change v{a) in an infinite sam- 
ple. It is obvious, however, that in order to "erase the 
memory" of short-range disorder, many more RG itera- 
tions have to be performed or, equivalently, much larger 
system sizes should be analyzed. Moreover, the larger 
the ratio Wo/Vb, the more challenging the numerics be- 
comes. At this point, we emphasize that in quantum 
percolation the short-range term emulated by the ran- 
domness in the phases has a huge magnitude. Indeed, if 
we choose in the first RG step all five SP's to be iden- 
tical with power transmission coefficients equal to 0.5, 
then due to the phases, the width of tiie Q{z) distribu- 
tions after the first step is already ±2.5.E3 This translates 
into an enormously wide spread in the transmission co- 
efficients of effective SP's ranging from 0.075 to 0.92. In 
order to suppress this intrinsic "quantum white noise" , 
one either has to perform more RG steps or to increase 
the magnitude of Vq (parameter (j in the notation of Sec. 
III). Both strategies are limited by numerical instabili- 
ties. In particular, a larger (3 leads to more weight in the 
tails of the Q{z) distribution in which the uncertainty 
is maximal. In other words, at large (3 the role of rare 
realizations is drastically emphasized. 

In fact, if we had to draw a quantitative conclusion on 
the basis of the accuracy we have achieved, we would base 
it on the curve in Fig. |^ corresponding to the maximal 
value /3 = 4. Actually, for this /3, the agreement with the 
extended Harris criterion is fairly good. In particular, for 
a = 0.5, we find « 3, whereas 2/a = 4. 

We also want to point out that the limited number 
(eight) of RG steps permitted by the numerics neverthe- 
less allows us to trace the evolution of the wave func- 
tions from microscopic scales (of the order of the mag- 
netic length) to macroscopic scales (of the order of 5^m) 
which are comparable to the sizes of the samples used in 
the experimental studies of scaling (see-a^., Refs. || and 
^ and much larger than the samplesOo used for the 
studies of mesoscopic fluctuations. 



C. Concluding remarks 

It was argued for a long time that the enhanced value 
of the critical exponent v extracted from the narrowing 
of the transition region with temperature has its origin in 
the long-ranged disorder present in GaAs-based samples. 
To our knowledge, the present work is the first attempt to 



quantify this argument. We indeed find that the random 
potential with a power-law correlator leads to the values 
of V exceeding v 2.35, which is firmly established for 
short-ranged disorder. Another qualitative conclusion of 
our study is that the spatial scale at which the exponent 
V assumes its "infinite-sample" value is much larger in 
the presence of the quenched disorder than in the un- 
correlated case. In fact this scale can be of the order of 
microns. This conclusion can also have serious experi- 
mental implications. That is, even if the sample size is 
much larger than this characteristic scale, this scale can 
still exceed the phase-breaking length, which would mask 
the true critical behavior at the QH transition. 

Our numerical results demonstrate that when the crit- 
ical exponent depends weakly on the sample size (large 
n in Fig. ^), the "saturated" value of v depends cru- 
cially on the "strength" (3 of the quenched disorder. Thus 
it is important to relate this strength to the observable 
quantities. We can roughly estimate (3 assuming that 
themicroscopic spatial scale (lattice constant) is the mag- 
netic length, while the microscopic energy scale (the SP 
height) is the width of the Landau level. We denote by 7 
a typical fluctuation of the filling factor within a region 
with size L. Then the estimate for /9 is /? ~ 7 (L/^)"^^. 
Naturally, for a given 7, the larger values of a correspond 
to the "stronger" quenched disorder parameter (3. 

Note, finally, that throughout this paper we have con- 
sidered the localization of a single electron. The role of 
electron-electron interactions in the scalin g of the integer 
QH effect was recently addressed in Refs. k3 and pi]. 
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FIG. 2. Top: P{G) (thin lines) as function of conductance 
G at a QH plateau-to-plateau transition. Symbols mark every 
20th data point for the initial distribution (■), the FP (□) and 
the distribution for RG step n — 16 (A). The vertical dashed 
line indicates the average of the FP distribution. Bottom: 
Corresponding plots for the distribution Q{z) of SP heights. 
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FIG. 3. Critical exponent obtained by the QH-RG ap- 
proach as function of effective linear system size L = 2" for 
RG step n. The error bars correspond to the error of linear 
fits to the data. The dashed line shows v = 2.39. Inset: v 
is determined by the dependence of the maximum Zmax.n of 
Qn[z) on a small initial shift zq. Dashed lines indicate the 
linear fits. 
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FIG. 5. Critical exponent u obtained by the QH-RG ap- 
proach as a function of RG scale 2" for /3 = 2 and different 
correlation exponents a. The dashed line indicates v — 2.39, 
which is the value that we obtain for uncorrelated disorder. 
For clarity, we show the errors only for a — 0.2 and 4. Inset: 
v vs 2°'"^^ does not scale for, e.g., 13 — 2. 
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FIG. 4. Moments {{G - (G))") of the FP distribution 
Pc(G) (o). Dashed lines are results from Ref. The dotted 
line corresponds to the moments of a constant distribution. 
Inset: Higher moments of Pc{G) following an exponential be- 
havior. The agreement of data and fit demonstrates the qual- 
ity of the fit result. 
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FIG. 6. Dependence of the critical exponent v on cor- 
relation exponent a for different /3 — 1,2,3, and 4 as ob- 
tained after eight QH-RG iterations. The dotted line indi- 
cates Qc = 0.84, following from the extended Harris criterion 
(Ref. M) for classical percolation. 
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